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ABSTRACT 

This work is a continuation of plasma oscillation studies of rare gas 
plasmas extending over the last five years at the U. S. Naval Postgraduate 
School. A low pressure Neon plasma was studied Electronic methods of 
obtaining continuous graphical plots of probe current squared, and first 
and second derivatives of probe current, with respect to voltage were used. 
Data derived from probe measurements is presented graphically in the form 
of: plasma density, space potential, average electron energy, and oscilla- 
tion intensities, as a function of probe position. Plasma densities derived 
from the positive ion section of probe current and a Druyvesteyn analysis 
are correlated with densities derived from detected electron oscillations. 

ae 
A correction factor on the order of le is shown to be necessary to 
Te 

correlate the density derived by positive ion current with that found by 
other methods. Electron oscillations were detected on all data runs through- 
out the length of the discharge, unlike the localized oscillations found in 
the argon plasma investigated the previous year. Ion oscillations corre- 
sponding to those in argon were detected throughout the length of the dis- 
charge. The detected electron oscillations corresponded to the Langmuir 
plasma frequency. Detector current from these oscillations increased by 
a factor of as much as 107 on decreasing plasma density gradients indicat- 
ing large growth in amplitude and decrease in frequency. The detected ion 
oscillations corresponded to the predicted Langmuir electrostatic sound 
waves and were found to peak in amplitude with the electron oscillations. 
These ion oscillations were assumed to have a wavelength nominally 


li 


the length of the plasma column (9cm), and a velocity equal to Yk Te ; 
177 ¢ 


These values of velocity and wavelength yielded frequencies of 94 KC 
and 75 KC in reasonable agreement with the 100 KC and 67 KC detected in 
a 75 ma and 65 ma discharge respectively. A tentative relation between 
the effects of ion mass and ionization potential on oscillations and dis- 
charge ionization phenomena is made suggesting that plasmas of light 
mass elements with atoms of high ionization potential may develop very 


strong electron and ion oscillations. 
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1. Imgroduction. 

In the years since World War II, the field of plasma physics has 
grown rapidly and with ever increasing importance. The study of plasmas 
has brought new knowledge into the fields of thermo-nuclear processes, 
radio communications, electronic circuitry, astrophysics, space propul- 
sion, ballistic re-entry phenomena, and the views of science on cosmo- 
logy. While these subjects are by no means the limit of its influence, 
they serve to illustrate the wide ranging effects of this specific field of 
basic research. 

The term plasma was originally defined by Langmuir and Tonks (4), 
and they stated: 

The word 'plasma' will be used to designate that 
portion of an arc type discharge in which the densities 
of ions and electrons are high but substantially equal. 
It embraces the whole space not occupied by the sheath. 

Plasma is not a rare manifestation of matter as might first be expect- 
ed by a casual glance at our immediate environment, but embraces very 
probably 99% or greater of the total matter in our universe. It is only in 
the dense matter environment such as laboratory facilities on earth where 
its existence must be so carefully sustained It exists naturally in the 
universe from the fringes of our ionized upper atmosphere throughout 
interstellar space, and the unbelievably energetic plasma phenomena oc- 
curring on the sun directly affect our lives on this planet. 

Plasmas generated under careful laboratory conditions have received 
increasing attention in the field of thermo-nuclear fusion devices. Initial 


attempts to generate a sustained thermo-nuclear reaction have met with 





repeated failures ina succession of experiments, both in this country 
and overseas. It is becoming increasingly apparent that an understand- 
ing of the simple basic behavior of plasma and its inherent instabilities 
is a necessity before more complex procedures of containment can be 
applied to it. 

It is toward this goal of better understanding of the basic phenomena 
in plasma physics that this particular research is carried out. This pro- 
ject is a continuation of the group work guided by Professor N. L. Oleson 
in the field of rare gas plasmas. In particular it is an investigation into 
the electron-ion oscillation instabilities that occur in a low pressure Neon 


discharge. 





2. Previous Work. 

Probe study of the characteristics of low pressure gaseous discharge 
tubes was placed on a firm theoretical foundation by Langmuir and Mott- 
Smith during 1924 in their analysis of a mercury vapor arc using probes 
(or collectors) of various geometries. (1) 

Further studies by Langmuir (2) in 1925 showed that three separate 
energy groups of electrons were apparently present in the discharge. He 


classified them as “primary electrons," “secondary electrons," and 


"ultimate electrons." The primary electrons were those which kept nearly 
all the momentum gained in the cathode ion sheaths. The secondary and 
ultimate electrons were two groups with Maxwelliam distributions of energy, 
the secondary's distribution proportional to the primary’s energy and the 
ultimate electron's distribution independent of the primary’s and lower than 
the secondary's. This indicated that some energy transfer mechanism 
existed in the discharge which Langmuir called scattering. Dittmer (3) 
disproved the hypothesis that a new mechanism, such as the possible large 
increase in the target area of excited molecules was the cause of scatter- 
ing, by showing that the scattering reached a constant value at some dist- 
ance from the cathode. Dittmer found that the region of scattering approach- 
es the cathode as the tube current is increased and predicted that high 
scattering of electrons might be due to oscillation in the plasma. However, 
he was unable to detect such oscillations with the equipment available at 
that time. 


Tonks and Langmuir (4) provided a simple theory for electron and ion 





=, : 
2] Ss > _ 
——> § aae= 
— me mma 
Om me ~ 
>> © 


> <a ¢ —_— 


— 
—_ = «& 
e—_— aE <a 
=. an 
en hs 
— 





oscillations, by treating the electrons as simple oscillators held ina 
jelly-like substance being the ions of the plasma. They predicted elec- 
trons oscillations having a natural frequency; 
ese 
We = (BE )* 
2 TC Wee 


and ion oscillations with an upper frequency limit of; 


} 
LW: = a ay 
Tentative verification of the theory was obtained by detecting oscillations 
which corresponded reasonably well with those predicted. 

Merril and Webb (5) in their investigation in a low pressure mercury 
arc discharge, using probe techniques, found scattering occurred in narrow 
regions, and was associated with oscillations, agreeing within ten per cent 
of the Langmuir frequency. Their conclusions were that scattering was due 
to oscillations receiving energy from the fast electron group. Continuing 
these investigations in the period 1947 to 1954, Armstrong, Emeleus, and 
Neill established that (a) oscillations are reproducible, (b) oscillations 
varied in frequency and magnitude over tube characteristics (6) and (c) the 
frequency of oscillation, the Langmuir frequency, was independent of probe 
bias over a large range, and varied directly with gas pressure (7, 8). 

Emeleus (9) demonstrated the value of visual observations ina typi- 
cal low voltage argon plasma in the micron pressure range. A sketch il- 


lustrating the major sectors follows: 








ow 


Fig. A. C, cathode; |. cathode sheath, 
2, primary beam; 3, meniscus, 
4, 5, scattered beams; A, anode 

Investigation of probe data during discharge studies by Easley (10) 
demonstrated the importance of probe cleanliness for accurate measure- 
ment technique. 

Oleson and Found (11) have shown that nonlinear log plots of probe 
data were obtained in discharges with striations present. 

It became increasingly evident that new methods of determining 
actual energy distributions in low pressure discharges was necessary 
Since the assumption of a strictly Maxwelliam energy distribution was 
not at all times valid. In 1939 Druyvesteyn (12) demonstrated that the 


electron energy distribution function, NgeY) could be found from the fol- 


lowing formula: 


Vp (eV) = — (ze)? dD * r 


Where N p(eV) is the number of particles of energy per unit volume, and 
in this work in units of reciprocal electron volts and reciprocal cubic 

: -] -3 . 
centimeters, i.e. eV cm ;A is probe area, and mg, e, Ip and Vp 
represent the quantities tabulated in the Table of Symbols and Abbrevia- 


tions. 








Waymouth (13) illustrated the utility of this expression for use in cases 
where departure from normally assumed energy distributions exists. 

R. H. Sloane and E. I. R. MacGregor (14) andR. L. F. Boyd and 
N. D. Twiddy (15) utilized modifications of a suggestion by Emeleus for 
obtaining the energy distribution. It was shown that a small a-c signal 
on the probe would yield a D. C. increase directly proportional to the 
second derivative of the probe current-voltage relationship. 

Boyd and Twiddy (15) also demonstrated that for a non-isotropic dis- 
charge a cylindrical probe will give an energy distribution which will attain 
negative ae between major electron groups. 

An R-F technique of determining plasma frequency was discovered by 
Takayama and Ihegumi (16) who found that if an R-F signal is superposed 
on a probe maintained negative to space potential that a marked resonant 
increase appears in probe circuit current at the Langmuir frequency 

Allis (17) and D. W. Mahaffey (18) have correlated plasma oscillation 
growth and decreasing plasma density gradients. Mahaffey reported an 
eighty-five per cent correlation between these parameters in a micron 
pressure mercury arc discharge. Putman and Collins (19) in previous work 
with N. L. Oleson also demonstrated the connection of these parameters 
in a low pressure argon discharge wherein peak oscillation intensity 
occurred along a decreasing density gradient. 

In 1949 Bohm and Gross (23) proposed a theory for plasma oscilla- 
tions which postulated that a longitudinal traveling potential field could 
be excited by a beam of high velocity electrons in transit through the 


plasma. This in turn would provide a mechanism for energy transfer from 





the beam to the plasma. 

Looney and Brown (24), however, in attempting to utilize this mech- 
anism by using a second beam of electrons were unable to find the pre- 
dicted oscillations, but rather postulated that they had found standing 
waves inthe plasma. They attributed these waves to interaction of the 
beam with large sheaths on the beam electrodes 

Demirkhanov, Gevorkov, and Popov (25) have shown that these oscil- 
lations can be produced as predicted, and have attributed the mechanism 
causing the effects found by Looney and Brown to secondary electrons 
ejected by the impinging beam on the anode. These secondary electrons 
oscillate in a potential well in front of the anode from which ejected 

G. D. Boyd, L. M. Field, and R. W. Gould (26) have also shown that 
this transit mechanism for electron plasma oscillations does occur and have 
succeeded in exciting them with a velocity modulated beam of electrons 
and have demonstrated that they grow in amplitude as they wave through 
the plasma. 

I. Alexeff and R. V. Neidigh (27) in recent work at Oak Ridge National 
Laboratories have found correlations between acoustic theory and kilo- 
cycle oscillations in a spherical glass discharge tube. The predicted 
Langmuir electrostatic ion sound waves are believed to have been success- 
fully observed with the following characteristics: 

a. Sinusoidal waveform. 
b Higher frequency modes of oscillation not integral 
multiples of the fundamental theory (This is in 


agreement with the boundary conditions of acoustic 





theory. 
| -l/2 
. Frequency proportional to Mi ; 
. Frequency proportional to p7! (D is tube diameter) 


. Frequency yields a value for electron temperature (Te) 


in agreement with values obtained by the other methods. 





3. Anomalies of Neon. 

In extension of the continuing project of studies of low pressure 
gaseous discharge phenomena, it was suggested by Professor N. L. 
Oleson, that a Neon plasma be investigated, to further extend data in 
this field. 

Prior to introducing Neon into this system Argon had been utilized 
for a test discharge and curves, as shown in Fig. 7, quite typical of 
those found by Putnam and Collins (19) for current-voltage characterist- 
ics were traced Slightly higher values of current were found due to the 
larger surface area of the cathode presently installed 

The system was thoroughly baked twice at 350 degrees Centigrade 
and flushed with Neon twice prior to taking any tube characteristic curves 
and making a data run. After the second bake out and flushing, the system 


7mm. of (Hg). The system was then filled with 


was evacuated to 2 7x10- 
9.6 microns (Hg) of Neon at room temperature, and the cathode brought up 
to operating temperature. Pressure then stabilized at 10 microns. 

It became immediately apparent upon taking the initial tube current- 
voltage characteristic curve that Neon was behaving considerably differ- 
ent in the same system than had Argon. Discharge currents were lower 
throughout the operating range for similar voltages, and severe instability 
occurred in the discharge at currents higher than 75 ma., at 10 microns 
pressure, and at lower currents at lower pressures. At currents below 30 
ma. at 10 microns pressure violent relaxation oscillations occurred, Fig.8, 


which sporadically would extinguish the discharge completelv. Voltage 


had to then be increased until re-ionization occurred. 





The relaxation oscillations were dampened as much as possible by 
placing a large resistance next to the cathode in the discharge current cir- 
cuit, anda large capacitor across the tube itself The values of these 
circuit elements were adjusted until an optimum damping effect was reach- 
ed. The values used for these components were 5500 ohms and 10 micro- 
farads. 

The instability occurring at the higher current range could not, however, 
be diminished except by going to considerably higher pressures. Ata 
pressure estimated at 50 microns - 70 microns it was possible to maintain 
a stable discharge up to 150 ma., however, since no accurate method was 
available of measuring this pressure, so that a return to the exact pressure 
could always be assured, it was decided to work in the more limited 10 
micron range. 

This high current discharge instability occurred in the following manner, 
and its threshold was always reproducible at ~ 75 ma. at 10 micron pres- 
sure. When 75 ma _ was reached, the voltage would slowly climb from its 
initial value to about 60 volts, then it would jump almost instantaneously 
to 325 volts. The current would decrease ina similar manner, slowly at 
first then almost instantaneously to a value of 15 ma., a value below the 
thermionic emission current of the cathode in high vacuum (22 to 35 ma.). 
The discharge in this form appeared as a narrow beam in the center of the 
tube with an apparent focus, independent of anode-cathode spacing, at 
about 40mm. from the cathode. Fig. 9. 

It was first thought that this might be due to the focusing effect of the 
slight cathode face curvature, however, when the anode was placed within 


10 





10 mm. of the cathode or less, the phenomena still occurred, so this 
hypothesis does not seem to be borne out. Further investigation of the 
literature revealed that Neon has been used extensively in slow sweep 
cathode ray tubes for “gas focusing" and is most effective in the 4 
microns - 6 microns range, Cobine (20). The 4 microns - 6 microns 

range is where the most severe form of this instability occurred, which 
we termeda “collapsed discharge," and tube currents of 40 ma. to 50 ma. 
were the maximum which could be sustained without collapsing. It is 
therefore believed that this effect was the primary reason for this upper 
limit instability in the Neon plasma. 

A cursory investigation by oscilloscope in this discharge form was 
remarkable for the absence of any oscillations and extremely low noise. 
The cavity wavemeter was also unable to detect any high frequency oscil- 
lations. Visual observations showed only a minor amount of the typical 
red glow associated with the normal Neon discharge illustrating low excita- 
tion and ionization occurring in this form. 

Further search into the literature while investigating this phenomenon 
revealed some other anomalies of Neon which it is felt deserve attention. 

It has been mentioned previously that lower currents were obtained 
throughout the current-voltage characteristic curves for Neon when com- 
pared with Argon. This at first seemed easily attributable to the higher 
ionization potential of Neon, 21.5 V, vice that of Argon, 15.7 V, however, 
it further became apparent that it was remarkable that such a comparable 


discharge did in fact occur at all. 
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Neon exhibits little Ramsauer-Townsend effect throughout a voltage 
range of 0 to 100 volts, and when compared with Argon its electron eo:li- 
sion probability (P;) is much smaller in our system's operating range. 
Fig. B. from R. B. Brode (21) illustrates this effectively. 

More explicitly, its relative ionization in comparison with Argon can 


be seen to be quite small as illustrated in Fig. C. from Von Engel (22). 
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and although this data was taken at 1mm pressure, we believe an order of 
magnitude extrapolation to a lower pressure is consistent for comparison. 
Also, the mean free path of an electron in Neon at a pressure of 10 mi- 
crons is, from kinetic theory, 9. 76 cm. a value which should be reasonably 
accurate since Neon having a tightly bound closed outer electron shell ap- 


proaches the classic spherical atom. This would lead one to believe that 


12 


the ionization of Neon would be considerably dependent on the cathode- 
anode spacing of the discharge up to 5 cm. or so, and relatively inde- 
pendent at greater distances. 

However, none of these effects were apparent in our discharge tube 
in the range expected. The currents as stated before were slightly lower 
than for Argon, but only by a factor of a few per cent at the most. The 
current and voltage were completely independent of cathode-anode spacing, 
as near as one could tell by reading a well calibrated milli-ammeter. for 
distances of 10 mm. to 100 mm. 

We therefore conclude that ionization is occurring rapidly within a few 
millimeters of the cathode and that the plasma may exist because of its 
oscillations interaction with the electron beam immediately upon leaving 
the cathode. As discussed later in the conclusion of our investigation. 
oscillations were found up to at least 2 mm. of the cathode face, which 


would certainly tend to confirm this hypothesis. 
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4. Equipment. 

Figs. 1 through 3 are photographs of the laboratory equipment arrange- 
ment. The equipment used in the experiment can be classed into four 
general headings: 

1. Gaseous discharge tube 

2. Vacuum system 

3. Frequency detection and measuring equipment 

4. D. C. electrical circuits. 
Discharge Tube. The tube envelope was locally manufactured using Pyrex 
glass and salvaged tungsten-uranium glass seals. Tube design criterion 
was based upon having a sealed system which could be subjected to bake 
out and, at the same time, allow a maximum amount of variation in the 
geometrical parameters 

The electrodes were all mounted on cylindrical iron cores sealed in 
pyrex to enable the electrodes to be moved within the tube by application 
of an external magnetic field. 

Electrical leads to the electrodes were coiled multi-strand oxygen-free 
copper, insulated where necessary with ceramic beads. Attachment to the 
tungsten seals at both the tube and electrode seal was by means of spot- 
welded nickel straps. 

The cathode was an indirectly heated dispenser type cathode supplied 
by Semicon of California, Incorporated. This cathode has an evaporation 
rate of one microgram per centimeter per hour and negligible sputtering. 
This gives the dispenser cathode considerable advantage over the conven- 


tional oxide coated cathode as regards tube contamination in this type of 
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work The cathode was 28 5 millimeters long, 22.1 millimeters in dia- 
meter, and had a concave face with a radius of curvature of 36.7 milli- 
meters. The heaters for the cathode were locally manufactured from 20 
mil malleable tungsten wire wound into a coil of approximately 15 turns 
and 14.2 millimeters in diameter. To prevent shorting between adjacent 
turns, the cathode was lined with two close-fitting alumina rings and an 
alumina disk at the face (see Fig. 5). To obtain the operating tempera- 
ture of 1075°C, required approximately 12.5 amperes. About one-half 
hour was needed to reach this temperature equilibrium. 

The cathode was mounted to the transport by means of a molybdenum 
jacket spot welded to four tungsten support rods (see Fig. 5). The rear of 
the cathode was closed witn molybdenum to minimize darkening of the tube 
by depositions of alumina and tungsten. 

The anode was a 26 millimeter diameter, 20 mil thick nickel disk sup- 
ported by a tungsten rod enclosed in a pyrex sheath. 

The probes were constructed of 5 mil tungsten wire sheathed in pyrex 
to within 3.9 millimeters of the end. To prevent conduction on the surface 
of the glass due to contamination, the probe has a 5 mil stand off from the 
glass sheath (see Fig. 5). The double probe was not used due to breakage 
in the final stage of assembly. 

Construction Problems. Early in the work, it was necessary to open the 
discharge tube due to an open cathode heater. This operation started a long 
chain of difficulties which terminated in a virtual re-building of the tube. 

The original heaters were constructed of 10 turns of 15 mil tungsten 
with an overall diameter of 17 millimeters. These heaters were coated with 
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a thin alumina coating to prevent shorting between adjacent turns. Several 
of these heaters were used; however, it was found that their useful life was 
quite short. It appeared that at the temperatures required for operation, the 
alumina would flake off and adjacent turns would short as the heater sagged. 
The next attempt to solve the heater problem was to use an alumina support 
on the center leg of the heater to give some additional rigidity to the heater. 
Somewhat greater success was obtained with this configuration While some 
shorting did occur, it did not appear to be progressive. Semicon of Californ- 
ia donated the present cathode, which in its original configuration had its 
heater potted in cast alumina. This would seem to be an ideal configuration, 
but the structure appears to be extremely fragile, perhaps too fragile, for an 
experimental arrangement. After two failures with this configuration and 
upon consultation with Semicon, the configuration described in Fig. 5 was 
used. 

In early attempts it was found that only very low emission currents 
could be obtained from the cathode. It was then discovered that the cath- 
ode was extremely sensitive to chloride poisoning, contamination by low 
vapor pressure metals, and contamination by water vapor and oxygen. 

The majority of these contamination problems were solved by proper 
activation of the cathode. The chloride poisoning, however, was irre- 
versible, destroying the usefulness of the cathode in general; if localized, 
it disturbed the emission in the vicinity of the poisoning. 

Activation of the cathode should be attempted after bake out. The pro- 


cedure is basically to bring the cathode to operating temperature (1 1602C) 
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in a high vacuum ~1075 millimeters of Hg, and apply a potential from 
anode to cathode of approximately 500V. As the cathode becomes acti- 
vated, a tube discharge current will be observed which will increase with 
time until activation is completed, at which time it will become constant. 
This procedure will require about two hours, depending upon the degree of 
cathode contamination. 

As mentioned previously, chloride contamination results in permanent 
damage to the cathode. The migration of other metallic ions to the emitting 
surface may also produce permanent damage to the cathode. It is there- 
fore necessary to keep low vapor pressure metals remote from the hot cath- 
ode. The materials recommended for cathode mounting are tungsten, tanta- 
lum, or molybdenum. Early attempts were made using a tantalum anode 
with the cathode heat shield and mounting of tantalum. On completion of 
bake out, it was noted that the tantalum anode was plated with a rather 
heavy coating of cupric, cuprous oxides and copper. The other tantalum 
surfaces were plated with a thin copper coating, in addition a powdery 
deposit of copper was found on the glass envelope (see Fig. 45, a,b,c). 
The only possible source for the copper was from the multi-strand copper 
wire. In view of these observations, it was postulated that the tantalum 
acts as a copper "getter" lowering the partial pressure of copper in the 
system, thereby permitting additional boil off of copper. If it were ab- 
solutely necessary to have copper and tantalum simultaneously in the 
system, this problem could possibly be minimized by an extremely slow 
bake out, always maintaining a vacuum in the order of 107° millimeters 
Org . 
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The final configuration utilized a nickel anode with all cathode mount- 
ings of molybdenum. A heat shield was not mounted. Mechanical joints 
were utilized wherever possible due to difficulty in welding molybdenum 
and tungsten even when using platinum as a "wetting" agent 

All copper leads were spot welded and supported with nickel straps 
This procedure minimized the amount of heat applied to the copper wire and 
reduced the possibility of contamination from either silver solder or flux. 
Vacuum System. The vacuum system was capable of a vacuum in the order 
of 107? millimeters of Hg and tested leak free with the Veeco MS-9 Helium 
leak detector which was capable of detecting leaks as small as 1x 190718 
cubic centimeters per second. Pumping was done with a three stage air- 
cooled oil diffusion pump and mechanical fore-pump. The vacuum system 
was divided into two sections: "internal" discharge tube, and vacuum 
gauges; and the "external" oil manometer and gas bottles for neon and argon 
gas. The "internal" system was isolated from both the pumps and the 
"external" system by liquid nitrogen traps. All stopcocks were vacuum 
type lubricated by Type N Apiezon grease. The oil in both the manometer 
and diffusion pump was Octoil-S. The "internal" system was baked out 
by a portable oven, while the "external" system tis baked out with heat- 
ing tapes. The anode was out gassed by means of an induction heater. 

Vacuum measurement in the high vacuum region ome to 1077 milli- 
meters Hg was by means of a Veeco RG-75 ionization gauge. To mea- 
sure micron range pressures, a modification was made to the Veeco Vacuum 
Gauge RG 31A by installing an additional range resistor to the circuit to 


provide a 1073 millimeter scale. This additional resistor in the pressure 
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multiplier circuit extended the scale of the meter by serving as a meter 
shunt for the ion current from the RG 75 ionization gauge. This modifica- 
tion gave continuity when compared to the i074 millimeter scale. On the 
basis of the manufacturer's claim for the RG 75 tube of linearity to oe. 
we assumed linearity for the 1-10 micron range. 

A Westinghouse type 7903 high pressure ionization gauge was install- 
ed in the system, however, control circuit difficulties precluded its use 
during our work. This instrument operates in the range of from 1 x 1079 to 
5x 107! millimeters of Hg and offers considerable utility for further work. 
High Frequency Measurement. A cavity wave meter, Test set, Model AN/ 
UPM 2Z, frequency range 80-1220 megacycles, was used to detect electron 
oscillation frequencies by capacitive coupling with a wire probe. The dc- 
signal from the detector crystal of the cavity wave meter was fed toa 
Hewlett-Packard Model 425 A dc micro-volt-ammeter. Currents as small as 
19712 amperes could be detected with this instrument. 

In previous work concern had been expressed that at the frequencies of 
interest (k-mcs) range excessive attenuation of signal might occur due to 
the coiled leads from the probe and their path through the iron case of the 
transport. To investigate this problem, the probe assembly was removed 
from the tube and measurements taken to determine the losses’ A signal 
was fed to the probe from a Hewlett-Packard 612 A frequency generator 
(450-1230 mcs) losses were measured through a detector crystal anda 
Hewlett-Packard 415 Standing Wave Ratio Meter for various configuration 


of the single probe the following results were noted. 
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Fig. D. Probe Attenuation vs Frequency. 

The characteristics for the double probe was far more complex and 
indicated quite high attenuation. 

Ion oscillation were detected and measured with a Tektronix Type 543 
oscilloscope. 

The circuits for obtaining frequency data were simply coaxial lines 
from the probe to the wave meter in the case of electron oscillations, and 
to the oscilloscope through a 1000 cycle high pass filter for ion oscilla- 
tions. The high pass filter was necessary due to a high background of 
electrical noise in the laboratory. 

Probe Circuits. Two sources of power were avilable to the probe The 


time varying signal was supplied by a Hewlett Packard Low Frequency 


Function Generator, Model 202 A, which was set for triangle wave output 
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varying linearly with time. The output signal was set to 15 volts peak to 
peak excursion and fed to a gain ten amplifier in the Donner Analog Com- 
puter, Model 30! with problem board Model 30-6. The amplified signal was 
then biased by means of the initial condition power supplies of the comput- 
er. This signal was then supplied to the probe to obtain the desired positive 
and negative voltage excursions with respect to the anode. For a constant 
voltage power supply to the probe, two of the initial condition power sup- 
plies of the computer were connected in series and connected to the probe 
This method was used to outgas the probe by ohmic heating prior to taking 
each probe characteristic. The voltage outputs of the initial condition am- 
plifiers were calibrated with a Donner Model 50 Comparator Bridge. 

Probe current as a function of probe potential was passed through a 
resistance on the order of 500 ohms where a voltage proportional to probe 
current was obtained. (The effect on probe potential due to voltage drop in 
the resistance is less than one per cent and considered negligible.) This 
voltage was fed through two gain ten amplifiers on the computers. The 
amplified signal was nen supplied to the F. W. Mosely Autograph X-Y 


2 2] 
recorder Model 2-S' either: direct as Ip. =. or 1% vs Vp. 


the Donner Analog Computer, Model 30, has 10 operational ampli- 
fiers which operate in the range -100 to +100 volts with a long time stab- 
ility better than 0.5% of full scale value, maximum drift characteristic of 
4 millivolts per hour with a short term random drift less than + 2 millivolts. 
Five initial condition regulated d.c. power supplier, 0 to 100 volts are 
available on the computer. 


2The F.W. Mosely Autograph X-Y Recorder, Model 2-S which has a 
range on the horizontal axis of from 0-7.5 millivolts to 0-150 volts full 
scale. The vertical axis has a range of 0-5 millivolts to 0-100 volts full 
scale. Input impedences on both scales are in the order of 2 megohms. This 
equipment is accurate to 0.25% of full scale on all ranges when calibrated. 
Calibration drift is negligible. 
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The X-Y recorder was operating with a sweep period of 20 seconds. The 
derivatives were obtained by two cascade R-C differentiator circuits 

with a time constant of .025 seconds. The Ip was Obtained by supplying 
the amplified voltage signal to both the X and Y inputs of a Donner 
Electronic Function Multiplier, Model 33. This equipment would accept 
inputs of + 100 volts in either or both X and Y inputs and produces an 
output of -0.01 XY. Use of this equipment was discontinued on the second 


run due to alignment difficulties at the last moment before the run. 
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5. Procedures. 
General Procedure. Prior to making a data run, the system was taken to 
as low a vacuum as possible, in the order of 107° millimeters of Hg, and 
Bei baked out with the portable oven at 350°C fora period of 10 to 15 
hours. During this time, external tubing was heated to 350°C with heat- 
ing tape. Upon completion of the bake out, vacuums in the range of 
2X mos? millimeters of Hg were obtainable. The next step was to bring 
the cathode to temperature with the pumps still on the system Since it 
was found that the alumina in the cathode would absorb considerable 
amounts of gas, the following steps were taken: when the vacuum stabil- 
ized, the cathode heater was secured; after the cathode had cooled, a 
charge of neon was admitted to the system with the pump secured; the 
cathode was again brought to temperature and a discharge fired for from 
15 minutes to one half hour, at which time the system was again pumped 
down to about 10~’millimeters of Hg. This procedure of flushing and fill- 
ing would be repeated two or three times to minimize the amount of foreign 
gas in the tube. The tube was again filled with neon and allowed to pump 
down to 9.2 microns as measured with the Veeco ionization gauge. While 
the absolute error of the Veeco gauge modification described on page 19 
may be large, it is considered that the pressures are reproducible to with- 
in two-tenths of a micron. 

Upon completion of the filling process, the cathode was brought to 
operating temperature and allowed to reach equilibrium. One-half to three- 
quarters of an hour were allowed for equilibrium to be achieved. 
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During this time, the electrodes were positioned with the catheto- 
meter using the following procedure: the cathetometer was set at its refer- 
ence value for the cathode. The cathode was then moved with the solenoid 
until its face was lined up with the reference mark of the cathetometer The 
cathetometer was then traversed to the probe position and this distance re- 
corded. Next the cathetometer was traversed to the desired anode position 
and the anode moved to coincide with the reference mark of the cathetometer. 
All measurements were made in the same direction to minimize backlash 
errors. Thee procedures should allow electrode placement accurate to with- 
in 0.1 millimeters. These same procedures were carried out, as applicable, 
for all subsequent electrode positioning. 

A discharge was then obtained with the oscilloscope across the tube. 
The value of resistance in the tube power supply circuit and capacitance 
across the tube were varied to minimize relaxation oscillations as observed 
with the oscilloscope. The X-Y recorder was then put in the discharge cir- 
cuit to record the tube characteristic, Vy asa function of Ig. With this 
information recorded, the X-Y recorder was placed back in the probe circuit- 
By 

The procedures from this point on became repetitive for each of the 
desired probe positions and were as follows: 

A. To observe ion oscillations, the anode was grounded, the os- 
cilloscope lead was directly coupled to the probe and the shield grounded 
common with the anode. Inthe first data run, the frequency and ampli- 
tude of the ion oscillations were read directly from the scope display. Dur- 


ing the second data run, all scope displays were photographed and analyzed 
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at a later time. 

B. Next the cavity wave meter was capacitively coupled to the 
probe and a frequency search conducted. The deflection of the micro- 
micro ammeter at resonance from its ambient value was recorded as a mea- 
sure of the electron oscillation amplitude. In many cases the deflections 
are very slight, a few per cent of the ambient current. It was therefore im- 
perative that the search be conducted very slowly to prevent resonance dips 
from being masked by random meter fluctuations which were not reproducible. 

C. The next step involved allowing the discharge circuit to float. 
At the same time, the probe circuitry which had been isolated from the dis- 
charge tube during frequency measurement, was then directly coupled to the 
probe. The output of the function generator was set for a triangle wave and 
adjusted to 15 voits peak to peak and fed to the X-axis of the X-Y recorder. 
The same signal was supplied to a nominal gain ten amplifier which was then 
calibrated to give a peak to peak excursion of 150 volts. This output was 
then biased by an initial condition power supply to, in the first run +43 to 
-107 volts and in the second run from +30 to -120 volts as measured be- 
tween probe and anode. (See Fig 5) 

The current drawn by the probe was passed through a resist- 
ance Rp to provide a voltage signal proportional to the probe current. This 
voltage signal was supplied to two gain ten amplifiers of the computer. The 
resistance Rp was adjusted to provide an output from the amplifiers of less 
than 100 volts, the maximum scale value for the Y axis of the X-Y recorder. 


The amplified voltage signal proportional to probe current was now plotted 
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as a function of probe potential. This same voltage signal was also fed 
through a pair of R-C differentiators 1n cascade. The outputs of this cir- 
cuit were then plotted as ae and 1p as a function of Vp. (See Fig. 6) 

Prior to each data plot, the probe circuits were isolated from the 
probe and a voltage of approximately +125 volts was applied to the probe to 
outgas it by ohmic heating. 

To allow for the differences in amplitude of the various signals, 
the calibrated scales of the X-Y recorder were changed, Ip vs Vp was plotted 
on 100 volts full scale and 10 volts full scale. The first derivative was 
plotted on 10 volts full scale, while the second derivative was plotted on 
I volt full scale. Attempts were made to use the 100 millivolt scale to en- 
large the plot for the second derivative, but the response was quite noisy 
and did not appear to be reliable enough for use. 

D. During the collection of probe data, the cathode temperature 
and heater current were monitored to keep the cathode within 5 °C of the 
selected operating temperature. Discharge voltage and current were record- 
ed for each probe position. During the first run, discharge current was 
continuously recorded in order to determine the effect on discharge current 
of drawing probe current. No discernable effect was noted even for maxi- 
mum values of probe current and its use was discontinued. 

E. At the completion of data collection for each probe position, 
the electrodes were repositioned and the balancing of the amplifiers check- 
ed. The entire procedure for a single probe position required approximately 


twenty minutes. 
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F. At the completion of each run, it was necessary to calibrate the 
data curves. This was done by computing the current necessary to produce 
the given voltage across resistance Rp. The first derivative curves were 
calibrated by measuring the slope of Ip vs Vp curve at a voltage correspond- 
ing to a maximum of the first derivative. The same procedure was carried out 
on the first derivative to calibrate the second derivative. The calibrations 
thus obtained were then checked against other sections of the curves, a 
good correlation was obtained. These calibrations are considered to be ac- 
curate at the maximum values to within + 5% of the true value as represent- 
ed on the plots. 

Data Run I. The first data run consisted of an investigation of a 75 milli- 
ampere discharge in neon by taking data every two millimeters over a length 


2 


of 40 millimeters. The data recorded on this run was Ip vs Vp, Ip*% vs Vp, 


a vs VD, ae vs Vp, electron oscillation frequency and amplitude, ion 
oscillation frequency and amplitude for each of the 19 probe positions. The 
purpose of this investigation was to determine the character of a neon dis- 
charge and to determine regions of maximum interest for further investigation. 
Data Run II. The second data run was taken ina 65 milliampere discharge 
and was set up primarily to determine the effect of a concave cathode face 
on the characteristics of the discharge, and to investigate the character of 

a discharge wherein the effects of Ramsauer free path are minimized, as in 
the case of neon. (See p. 12). 


The program consisted of taking probe data with the probe at 5, 10, 15, 


18, 21, 23, 28, 33, 35 millimeters from the cathode. For each of these 
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positions, the anode was to be located at 5, 20, 35, 50, 65 millimeters 


from the probe. The data recorded for this run was Ip vs Vp, Ip X 10 vs 


dip d2I — | 
Vp, Evo vs Vp, avez vs Vp, electron oscillation frequency and amplitude, 


ion Oscillation frequency and amplitude. 
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6. Results. 

Tube current-voltage characteristics curves are shown for data runs 
l and 2 in Figs. 10 and 1l. Both these runs were made at 10 microns 
pressure. Run 1 was conducted at 75 ma discharge current, andrun 2 at 
65 ma discharge current. A comparable curve for Argon is shown with that 
for Neon in Fig. 10. The various curves in Fig. 11 illustrate the absence 
of any effect due to increased anode-cathode separation from 10mm to 75mm. 
Only the general shape of these curves is reproducible since the power sup- 
ply voltage was changed by hand. The difference between them is attribut- 
able to this factor. It was found that if the voltage was stopped at any 
point, the current would reach an equilibrium value independent of Cathode- 
Anode spacing. 

The region of instability which eventually results in a collapsed dis- 
charge follows the rather short level section of the curve. The voltage rise 
in this part of the curve is clearly evident. A photograph of the relaxation 
oscillations occurring below 30 ma. is shown in Fig. 8., one of a collapsed 
discharge in Fig. 9. 

Typical probe characteristic curves for run 1 are shown in Figs 12, 13, 
14, and 15. They are for probe positions near the cathode (Fig. 12), at the 
beginning of the region of maximum interaction (Fig. 13), within the region 
of maximum interaction (Fig. 15). Probe characteristics curves for run 2 
with a fixed cathode-probe distance of 15.5mm and anode to probe distances 
of 5, 20, 35, 50, and 65mm are shown in Figs. 16, 17, 18, 19 and 20 re- 


spectively. These illustrate the gradual decrease then the increase in length 
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of the high energy tails on that section of the second derivative giving the 
energy distribution of the randomized electron group, and the gradual de- 
crease then the increase in size of that section representing the "beam” 
electrons. This result occurs as the anode-cathode distance is increased. 

The second derivatives of Ip with respect to Vp are een from these 
curves to approach zero more rapidly than they should. A changing slope 
may be seen on the curves representing the first derivative extending out 
to a considerably greater negative voltage than is indicated by the second 
derivative plot. Corrections made for this effect are discussed in our con- 
clusions. 

It was noticed from these curves, which were traced from a positive to 
negative probe voltage, that the inertia of the plotting arm of the X-Y re- 
corder causes a shift to the right, or toward to more negative voltage, for 
both the first and second derivative plots. This is illustrated by comparing 
the location of the maximum value of slope ofthe Ip vs Vp curve with the maxi- 
mum value of the first derivative, and the maximum value of the first deriva- 
tive with the zero of the second derivative. A shift along the horizontal 
scale of minus lvolt in each case is seen. 

The character of the representative curve in the region of maximum in- 
teraction (Fig. 13) illustrated the effect of the strong oscillations occurring 
at this point on the first and second derivative plots. 

Representative plots of space potential vs distance along discharge are 
shown in Figs. 21, 22, 23, and 24. Fig. 21 for the first data run shows the 


close agreement obtained between the Zero of the d2] , the maximum of 
dVp 
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dip, andthe "first break" in the In Ip vs Vp plot for probe curves taken 
dVp 
in the most "Maxwellian" section of the discharge, on the anode side of 
the maximum interaction region. 

Several electron energy distributions are shown in Figs. 25. 26 aand 


V2 "Ip vs probe retard- 
dVp 


ing potential directly from the probe curves after correcting for space poten- 


b, and 27. These were obtained by plotting Vp 


tial. The presence of the beam electrons at 45 to 50 volts energy is clearly 
seen in Fig. 25, a position 6mm. from the cathode inrun 1. Fig. 26 a and 

b illustrates the “smeared out" appearance of the distribution in the maximum 
interaction region, 30 mm. from the cathode, and Fig. 27 shows the random- 
ized "near Maxwellian" distribution well past the interaction region 38 mm 
from the cathode. These latter two positions are also from run 1. 

There is an absence in all these distributions of any intermediate or 
very high energy electron groups. One distribution from run 2 is includedin 
Fig. 28. This was the only distribution which shows three distinct energy 
groups ina total of 64 probe positions in both runs. 

Electron and ion densities are shown in Figs 29 a andb, 30, 31 and 32, 
for run 1, and three groups of 40mm, 55mm and 70 mm anode-cathode spac- 
ing from run 2, as a function of distance along the discharge The plots 
show the densities predicted from frequency, the Druyvesteyn analysis, 


2 


and the slopes of the positive ion Ip“ vs Vp portion of the probe curve. The 


last method's results are shown in run 1 both corrected and uncorrected for 


the ratio of 10n to electron temperatures, ii. 
Te 


The correlation between electron density and oscillation (or detector 
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current) intensity 1s shown in these figures wherein the density tn this 

case is derived from the frequency detected. The maximum error in de- 
tected frequency is less than 2 megacycles in the general range of our 
operation with the AN/UPM2 cavity wave meter. The maximum error in de- 
tector intensity is +3 units of the various scales used. In general both of 
these are within the plotting circles on the graphical plots. A plot of den- 
sity (derived from frequency) with space potential (Vsp) and electron tempera- 
ture (Te) is shown in Figs. 33, 34, 35, and 36. 

Photographs of ion oscillations varying from 50 KC to 66.7 KC are shown 
in Figs. 37, 38, 39, 40, and 41 for run 2. A possible higher frequency com- 
ponent of 300K (is dec2rnable superposed on the lower frequency traces. No 
ion oscillations near the upper limit predicted by Langmuir (about 2.3mc for 
our discharges) were detected. Plots of the ion oscillation amplitude vs 
distance are shown together with plots of electron oscillation intensity and 
density, in Figs. 42, 43, 45, for both runs 1 and 2. 

Electron oscillations were detected throughout the discharge in all 
runs. The oscillation frequency in general varied from 700 megacycles 
near the cathode to a predominant 500 megacycle range throughout the rest 
of the tube. Usually a slightly higher frequency was found near the anode, 
and several other frequencies would appear sporadically on the anode side 
of the maximum interaction region. Possible explartation for these frequen- 
cies are covered in section 7 of this work. 


A tabulated result of all frequencies detected in run 2 follows 
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Cathode Probe to Electron Electron Ion Osc. Ion Osc. 
to probe Anode Osc.Freq. Osc.Amp. Freq.Fi Amp. Vi 


(mm.) (mm.) Fe (mc) Ic amperes’ (kc) (mv) 
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5 35 702 s5x1071! 62.5 26 
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Cathode Probeto Electron Electron lon ‘Ose « Ion Osc. 


to probe Anode Osc.Freq Osc.Amp. Freq. Fi Amp Vi 
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: : a 1x 107° : 2 
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Cathode Probe to Electron Electron 1OneOsc. Ion Osc. 


to probe Anode Osc.Freq. Osc.Amp. Freq. Fi Amp. Vi 
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7. Interpretation of Results. 
General 

It was found in preliminary work that analyzing space potential, elec- 
tron density, and temperatures by those methods of Langmuir which assume 
Maxwellian distribution of electron energies produced an extremely poor 
correlation with data obtained through the several other methods. In view 
of the general non-Maxwellian behavior of the electron energy distribution 
function, and the Ip vs Vp curves, it was decided to rely upon those me- 
thods of analysis which are not predicated upon a Maxwellian type distri- 
bution of electron energies. 

Before introducing the data from our probe investigations of the dis- 
charge, its reduction, and our conclusions, we believe that some amplifi- 
cation of our statement concerning the existence of oscillations close to the 
cathode and the bearing of these oscillations on the maintenance of the 
plasma in Neon is in order. 

We believe that these oscillations are connected with an early inter- 
action of the beam electrons with the gas after they leave the cathode 
sheath. Classical kinetic theory. and experimental results of electron 
ionization efficiency, discussed in section 3 of this paper, predict a much 
lower ionization of the Neon than that which we have found in this region, 
and also predict a definite increase in ionization with increased cathode- 
anode spacing out to a distance equal! to the mean free path. We found 
that this increase did not occur, and that a relatively constant value of 


plasma density existed throughout the discharge, differing by a factor of 
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three at the most. 

It is believed that a mechanism described as follows may account for 
this. 

The high energy electron beam emitted from the cathode will cause 
some ionization after leaving the cathode sheath. The beam will then con- 
tinue to interact with this rarified plasma adding kinetic energy to the elec- 
trons and ions through the mechanism described by Bohm and Gross (23). 
This would have two effects, first prevent recombination of ion electron 
pairs, and second, cause increased ionization by the plasma electron-gas 
atom collisions taking place inthe region. If true this would indicate that 
two complementing ionization mechanisms might occur in a gas, and help 
to explain the lack of oscillations found near the cathode by Putnam and 
Collins (19) in the Argon discharge. Argon would ionize far more readily 
due to the resonance effect (Raumsauer effect) near its ionization potential 
as seenin Figs BandC section 3, than does Neon. This rapid ioniza- 
tion due to electron collision, and Argon‘s increased mass over Neon, 
would mean that beam electrons traveling sufficiently close to an atom 
would cause more ionization by this collision effect and transfer less ener- 
gy to the ions due to their traveling potential field effect. Also, the Argon 
plasma electrons in this region may not be able to gain large amplitude os- 
cillations since they can more readily ionize their nearest neighbor atoms 
than in the corresponding case of Neon. 

Another factor entering into the absence of these oscillations in Argon 
could well be our extreme good fortune in capacitively coupling our wave- 
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meter to probe. We were able, after many tries,to find an extrerely sen- 
sitive coupling which allowed us to detect oscillation amplitudes down to 
190712 amps in detector current. These frequencies near the cathode were 
always very weak (ig7il amps) and a smal] change in coupling impedance 
could easily have prevented their detection. 

Probe Characteristic Curves. 

It is noted 1n inspection of the probe characteristic curves that the 
relative density of the energy distributions changes with increased anode 
cathode separation for a single probe position with respect to the cathode. 
In fact as seen in Fig. 48, the high energy beam group increases, aftera 
gradual decrease, with increased separation as does the randomized plas- 
ma electron group. 

We conclude from these effects that two mechanisms are occurring. 
The first which yields our initial high density with the anode close to the 
probe is the emission of secondary electrons caused by the beam imping- 
ing on the anode. Moving the anode away from the probe would cause 
the effect of these secondaries to decrease, since they are trapped in 
space within a few millimeters of the anode (25). 

The second mechanism which would cause the significant increase of 
these groups as the anode is moved farther from the cathode is not at all 
clear. We have several] hypothesis and they are: 

1. Possible convergence of peripheral beam electrons toward 


center of discharge as anode 1s moved further away 
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2. Possible influence of field lines of concave surface of cathode 
causing a slight focusing effect on beam when anode spacing 
is large. 

3. A high positive space charge found 20 to 30 mm from cathode 
in plasma may tend to attract peripheral beam electrons toward 
the center of the discharge. Fig. 21 illustrates an example of 
this charge and its location in the plasma. 

From the appearance of the probe curves in the region of maximum in- 
teraction as illustrated in Fig. 13, the difficulty of extracting information 
about electron energy distributions in this area becomes apparent. The 
effect of the beam plasma interaction on the second derivative of Ip with 
respect to Vp is so large, that no absolutely valid information could be 
taken by the Druyvesteyn method in this region. In all figures showing 
data plots where this effect occurred we have shown cross hatched areas 
It is noted that in general it is these regions where our data correlation is 
poorest. In other areas we believe the correlation to be quite good. The 
example probe characteristic curve shown in Fig. 13 is the most severe 
example of this disturbance. It was quite reproducible, and it can be seen 
that the second derivative was traced three separate times. The disturb- 
ance is also apparent on inspection of the Ip vs Vp ond gip vs Vp curves. 


Electron Energy Distributions. 


An analysis of the electron energy distribution was carried out using 


Life, 
the method of Druyvestyn (12). The product of Vp / d“Ip was plotted 


dV2 
p 
against Vp to obtain the shape of the energy distribution function. The area 
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under this curve could then be interpreted as the electron density by 


numerical integration where 


VBV) = (amet dy 


e 


where N = electron density in lectin ent. -eV 


A area of probe 


For our situation it was determined that 


N bev sae 2.7210 Vdd and NV = (wae) eV) 
r 


for all regions where the second derivative has positive finite value. The 
results of this procedure produced electron densities which were in error, f 
when compared with the densities obtained by the relationship N = (sa50) i 
by a factor of about 1.5 when compared at all points. Upon further inspec- 
tion of the probe characteristics curves, it became apparent that not all 
electron energies were being included in the energy distribution. By in- 
spection of the a vs Vp curves (Figs. 12 through 20), it can be seen that 
there is a recognizable slope which extends out to approximately -65 volts. 
At the same time, it can be seen that the effect of this slope does not in- 
dicate in the plot of the second derivative of current with respect to volt- 
age. We therefore concluded that perhaps a relatively large number of elec- 
trons which should have been included in our energy distribution had been 
missed. To obtain an approximation to the magnitude of this loss, we as- 


sumed that the trailing edge of the energy distribution (Fig. 27) from the 
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last readable value obtained from second derivative would continue out 

as a straight line to zero at -65 electron volts energy, rather than to where 
the second derivative had appeared to have gone to zero. This operation 
increased the value of the electron density by a factor of about 1.4. Having 
established why our electron densities were in error and realizing that this 
method of correction left a great deal to be desired due to the small magni- 
tude of the effects as plotted, we decided to apply a constant correction 

for what we felt was the mean distribution of missing electron energies. 


This correction was obtained by picking four densities obtained from 


2 2 
n -F) and comparing them to n= Kf? ve , where k is to be our 
gis 


empirically arrived at constant. This gave us a value fork of 4 4x1] 
and n was then considered to be n= 4.4 x 01 fyp/ 24g When this 
correction was applied to the remaining 37 energy distributions, a good cor- 
relation was obtained as noted in Figs. 29a, 30, 31, and 32. It can also 
be seen in these figures that the correlation is poorest in the region of 
maximum interaction where the second derivative of current with respect to 
voltage becomes difficult if not impossible to evaluate (see Figs. 13 and 
14). 

Plasma density was also obtained from an analysis of the positive ion 
region of the probe characteristics curves by the method of Langmuir (+). 


Ion densities are obtained from the relationship: 


nz. _ Wm: ASG) 
mes ZeA ol Vp 


4] 


isa - 
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ee 
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where nj; = ion density in ions/cm? 


i] 


A probe area in cm2 


mj; = mass of ion 
This analysis of the positive ion region does not assume a Maxwellian 
energy distribution as does the similar Langmuir method for finding electron 


é vs Vp curve in the electron region of the curve. 


density by the slope of the Ip 
The required slopes were obtained by squaring the values of Ip for every five 
volts Vp for the linear positive ion portion of the Ip vs Vp curves. When the 
ion densities were computed by this method, they were found to be about 
three times greater than the densities obtained by Druyvesteyn analysis and 
from the relationship n {t)° (Fig. 29a). L.S. Hall (28) has postulated 

a proportionality factor for the Langmuir relationship which is "expected to 
be in the order of a . We made attempts to bring our ion densities into 
better agreement with electron density by applying this proportionality fact- 
or “in the order of Tis . In computing the correction, we used for Te the 
mean energy of the randomized plasma electrons. For Ti, we assumeda 
constant ion temperature across the discharge of 1000°K, which seemed 


reasonable for our discharge. The next step was to correlate n (electrons) = 


kK Ti 

Af Ji: where for a best fit k = 3. This produced a relationship n 
e 

(corrected) = 


3h 
=—- i. 
We 


data points of run one and two, Figs. 29a, 30, 31, and 32 resulted. 


When this correction was applied to all of the 


It was considered significant that the agreement abtained between this 
method and the Druyvesteyn analysis agreed both in shape and magnitude, 
though showing some deviation from the densities obtained from n= soET 
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It was further considered significant that this correction removed the major 
fluctuation in ion density in the region of the anode. It was noted that the 
agreement was poorest in the region of maximum interaction as might be 
expected. 

To further investigate the proportionality factor, we analyzed the data 
of Putnam and Collins (19) which had also shown an ion density consider- 
ably higher than the electron density in Argon. Using the same technique 
and assumptions, we were able to obtain the same sort of correlation 
] 


described above by using a correction formula, n (corrected) = pa Nj- 
e 


Space Potential vs Distance. 


The plot of space potential vs distance for run 1 with a constant 
anode-cathode spacing of 40 mm shows an interesting effect which we cor- 
roborated by three separate methods. The data was taken from the zero of 
the second derivative of lp with respect to Vp, the maximum value of the 
first derivative of Ip with respect to Vp, and the first break from linearity 
of the In Ip vs Vp curve. The latter method, as stated previously, 1s not 
to be relied on too heavily because of our non-Maxwellian electron ener- 
gies, however, the gross effect of an increasing space potential some 
10 mm away from the anode is quite apparent In the most Maxwellian sec- 
tion of the discharge, i.e., on the anode side of the maximum interaction 


region, we find close agreement between all three methods of determining 


the space potentials. 


Pe, 
INote, it might be significant that yo = ay and that 


atomic weight of A is twice that of Ne e e 
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We have therefore decided to rely most heavily on the zero of the 
second derivative for our absolute value of space potential in all runs, 
and only this value was plotted for our 40mm, 55mm, and 70mm (anode 
cathode spacing) groups from run two. 

Unhappily in these runs (which were actually a collection of data from 
our investigation of the effects of various anode cathode spacings with a 
fixed cathode probe distance) we have few probe positions. However, we 
believe that the effect of increasing space potential some 30mm from the 
cathode may be indicated in them. There is a slight but definite increase 
in the second 40mm group and the 55mm group in Figs. 22 and 23. 

We are far from being positive about what causes this positive space 
charge effect, however, we postulate three mechanism which could ac- 
count for it. 

First, the cathode face curvature may result in some "bunching" of 
field lines in this area causing the ions to remain there, whereas the more 
mobile, more energetic electrons could escape. 

Second, the ion focusing effect of neon, mentioned in section 3. may 
be at a maximum in this area, causing the ions to collect in the region. 

Third, the beam electrons which have been randomized in direction in 
the maximum interaction region may be leaving this area toward the Tube 
walls as the “electron cloud" drifts down the tube toward the anode with 

Tinie method was brought to our attention by Professor K. G. Emeleus 
of Queens University, Belfast, Ireland, who very kindly allowed us to 


view an as yet unpublished thesis of A. Garscadden during a visit to our 
laboratory. 
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the normal electron drift velocity of the discharge 
Electron Oscillations. 

The correlation between electron oscillation amplitude (from detector 
current intensity) and electron densities as illustrated in Figs. 33, 34, 
35 and 36, we believe to be extremely significant. Fig. 33 illustrates an 
almost mirror image effect between these two parameters. Figs 34, 35 
and 36 also show remarkably good correlation. and we feel that these are 
of special value in corroborating this effect since they are taken froma 
collection of discharges with the same cathode-anode spacing, discharge 
current, and pressure, but over a period of several hours with two com- 
plete changes of gas in our system. 

This effect agrees with the theory that these oscillations should grow 
traveling down a decreasing plasma density gradient Putnam and Collins 
(19) also found this effect as did Mahaffey, et al (18). However, Mahaf- 
fey's group believed that many of the small intensity changes in this data 
might be due toa "frequency pulling effect" of the probe. Since Demuirk- 
hanov, Gevorkov, and Popov (25) were able to correlate many of their fre- 
quencies detected by probe, with signals detected by receiving equipment 
using an antennae outside their discharge tube, we believe that our fre- 
quencies detected through the probe are correct, and hence the densities 

f \4 
obtained from them through the relation n = (ssxa) are accurate. 

It should be noted that the intensity values are plotted on a rough 
logarithmic scale, and that small decreasing density changes yield very 
large changes in oscillation intensity. We conclude, therefore, that the 
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oscillations do start on peak plasma densities and grow rapidly down the 
decreasing density gradient, and that minor density changes have large 
effects on increasing and decreasing the amplitude of these oscillations 

We have also found, as illustrated in Figs 33, 34, 35, and 36, that 
the plasma density follows quite well the average electron temperature 
(plotted in electron volts) as calculated from the energies described by 
our distribution function. This result is to be expected, because the areas 
of higher electron energy should be able to maintain a more dense plasma 
The higher energy electrons will have less chance of 10n-electron recom- 
binations and aid in ionizing neutral atoms in this region. 

The result in turn gives us further evidence for the accuracy of our 
Druyvesteyn analysis of electron energy distributions, since we used this 
to predict our average energy, whereas the densities which are correlated 
with these energies come directly from the detected frequencies 

The method used to determine these energies was basically a numeri- 


cal integration of the electron distribution function in the following manner: 


iy = Fhe New sew 
JSNG@Y) dev) 


Electron oscillations were found, as tabulated in section 6 for run 2, 
to extend from cathode to anode in both data runs ] and 2. The increase 
in frequency found near the anode we attribute to the increased density of 
electrons resulting from secondary emission from the anode. Even very 


low energy electrons have a considerable secondary emission effect on 
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nickel (the anode material) 3, and the range of this effect agrees quite 
well with what would be necessary to produce the increased frequency 

in this region. The previously mentioned positive space charge on the 
anode side of our maximum interaction region would aid in removing these 
electrons from the anode upon their emission. 

Ion Oscillations. 

Ion oscillations were also found throughout the tube, and they were 
seen to peak in amplitude with the electron oscillation intensity. Photo- 
graphs of these oscillations show this effect (Figs. 37, 38, 39, 40, and 
41) as does the plot of oscillation intensity illustrated in Figs 42, 43, 
44, and 45 The amplitude of the ion oscillations, although peaking with 
the amplitude of the electron oscillations, shows none of the minor varia- 
tions which in the case of the electron oscillations are attributed to in- 
creasing and decreasing density changes. 

No ion oscillations were found in the 2.3 megacycle range as pre- 
dicted by the Langmuir relationship f, = fe (Me) 1/2 All oscillations 

i 
detected were of a much lower frequency, and we believe that these are 
the electrostatic sound waves predicted by Langmuir and Tonks (29) 
Using this assumption we investigated the boundary conditions on 


our plasma column in order to infer what mode of longitudinal oscillation 


might be expected. A visual inspection demonstrated that the most 


SH.S.W. Massey andE.H.S. Burhop, Electron and Jon Impact Pheno- 
mean, Clarendon Press (1952) Chapter V of this text gives a comprehensive 
discussion of this phenomena The curves for Copper having a value of 
maximum emission coefficient similar to Nickel, p. 306, illustrates the 
small but finite emission which may be expected from primaries of 40 volts 


and below. 
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probable situation was that of a free boundary surrounding the plasma on 
all sides. The glowing plasma column was surrounded by a dark space 
approximately one centimeter thick between itself and the tube walls. 
This space we believe consists primarily of neutral atoms within which 
ion-electron recombinations are occurring 

Simple acoustic analogy (30) indicates a longitudinal mode of vibra- 
tion for this boundary condition wherein the frequency is related to the 
length of the vibrating column by the expression: 


f= ne in =a) 


> 


] 
Zi Z| A 


where f is frequency, n an integer, c the wave velocity, and | the column 


length. Langmuir and Tonks predicted a wave velocity on the order of 
YKT. 


Mi 
tion the results are: 


es and when this expression is introduced into the frequency equa- 


21 Mi 


The above expression was used to predict frequencies from data runs 
1 and 2, and the following assumptions were mado. 

a. n= 2; since Langmuir indicated a low frequency cutoff wherein 
» was on the order of the discharge tube dimensions. 

b. 1= 9 cm: the approximate length of our plasma column. 

c. Y= 5/3; as used by Alexeff and Neidigh (27) 

cp iste: Boltzman's constant x electron temperature, the value used 
was actually the average value of (eV) converted to ergs at the 


location where ion oscillations were a maximum. For run 1, (eV) 
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= 9 electron volts, for run 2 (eV) = 6 electron volts. 
Utilizing these assumptions in the equation. with the appropriate 


value of the electron temperature for each run, gave the following results. 


= 
-/2 
ge = [yk Te f/f ZS “94/0 ergs = &.76 aioe om /see 
47}. 3.35 x /0-** ams 
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= 6:9) 70 am sa 
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Fron = = LL Le cm/sec — Se kilocycles 


We feel that these values agree reasonably well with the predominant fre- 


quencies detected in run l of ~ 100 KC, and those detected in run 2 of ~ 


See 6 79KC.. 
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8. Comelusions. 

In general, it is our belief that the following conclusions may be 
drawn from our work. 

First, the mechanism of oscillation generation does occur quite close 
to the cathode in our type of discharge tube, and may be dependent upon 
several factors. These factors are; the mass of the gas ion, the ioniza- 
tion probability of the atom, the ionization potential of the atom, and the 
plasma density in this region. 

In neon, wherever we placed the probe, we detected oscillations 
which correlated with the predicted Langmuir frequency of f = 8980 nlf2 : 
In Argon and other gases, however, these oscillations have been consider- 
ably more localized, usually in the vicinity of the meniscus or maximum 
beam-plasma interaction region. 

We postulate that the lighter mass of the neon atom, its low ioniza- 
tion probability, and its higher ionization potential allows the plasma 
electrons and ions near the cathode to oscillate with greater amplitude 
than in the case of such gases as Argon, where, in this region, the elec- 
trons, in gaining a few volts of energy, would soon cause other ioniza- 
tions rather than continue to oscillate and increase in amplitude of os- 
cillation. Thus, although oscillations would occur in this region, they 
wouldtend not to grow in amplitude, and be extremely weak in intensity. 


Even in Neon, our weakest signal was in this region and was of the order 


of io amperes. 
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It would be expected that once the beam electrons became random- 
ized with the plasma in the maximum interaction region that these effects 
would be minimized allowing detectable oscillations in all gases, since 
all the plasma electrons have a much greater average energy and would 
remain in an oscillating state in the ion matrix. 

The mechanism of interaction itself we believe to be that predicted 
by Bohm and Gross, wherein a traveling potential field is set up by the 
beam in transit through the plasma. 

Second, at the face or cathode side of the maximum interaction region, 
we see in our discharge complete randomization of the electron beam with 
the plasma. Wecall it the "brick wall effect" since it appears as if the 
beam electrons are suddenly stopped in their transit of the tube and scatter- 
ed isotropically. The intensity of oscillations in this region and further on 


in the interaction region reached values of lie amperes of detector cur- 


7 gil 


fentrathermthan*thes 10 “and 1 amperes elsewhere in the discharge. 
We postulate that these violent oscillations and their associated micro- 
fields are literally breaking up the directed beam. 

Third, we believe we have shown that plasma densities derived from 
slopes of the ion portion of the curve do indeed need to be corrected hy 
some function of The gross correction we made to our own data and 
that of Putnam and Collins’ work in Argon illustrated that reasorably good 
correlation with densities predicted from other methods could be obtained. 
Also, in the region where ion oscillations were most intense, and 10n 


energies greatest, use of our assumed correction factor based on our 
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average ion temperature, showed the poorest correlation with other densi- 
ties determined. 

Fourth, the 10n oscillations which we detected and photographed, we 
conclude are the electrostatic sound waves predicted by Langmuir. They 
agree well with his predictions and appear to be related by wave length 
to the longitudinal dimensions of our plasma column 

Fifth, we have found that in neon the anode-cathode spacing has 
little effect on frequency once the anode is sufficiently far from the cath- 
ode that beam electrons do not emit copious quantities of secondary 
electrons We do not know why the high energy beam electron group 
seems to grow with very large spacing but we have postulated three hypo- 
thesis in the preceding section. 

Sixth, we believe our high current instability effect as discussed in 
section three is due to ion focusing of the electron beam. It is an extreme- 
ly quiescent discharge with no detectable ion or electron oscillations occur- 
ring. We have recently put Argon in our tube and were able to draw up to 
700 milliamperes without the discharge collapsing, so we now believe 
even more firmly that this effect is an anomaly of neon. 

Seventh, and possibly the most important, we believe we have definite- 
ly confirmed the violent growth in amplitude of the electron oscillations on 
a decreasing plasma density gradient. Even the fine structure of the curves 
show this effect, and we feel it is predominantly real rather than a probe- 


plasma interaction. 
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9. Recommendations. 

In keeping with our observations and results, the following areas are 
recommended for further investigation: 

a. Continue this type of investigation in other gases in an attempt to 
determine whether ion mass. and ionization potential are major factors in 
establishing the intensity of electron and ion oscillations and their loca- 
tion in the discharge. An investigation of discharges at similar currents, 
pressures, and anode cathode separation using He, Ne, A, Kr, Xe, and 
mixtures of these should prove definitive in this investigation 

b. As an adjunct of the investigation recommended above, additional 
information should be obtained regarding the possible ion mass dependence 
on the form of the Ti/Te correction to the ion densities obtained by analy- 
sis of the positive ion region of the probe characteristics curve. 

c. Further investigation should be made of the low frequencies found 
in the kilocycle range to determine if they are longitudinal ion oscilla- 
tions as supposed, radial modes of vibration, or a coupled oscillation 
which might be a function of both sides. 

d. Additional techniques which are recommended include: use of a 
high sensitivity UHF receiver to be used in an attempt to detect oscilla- 
tions external to the discharge, and addition of a spectrographic analysis 


capability which will be extremely important, if mixed gases are used. 
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Fig. 8. Oscilloscope Photograph of 
typical relaxation oscillations 





Fig. 9. Photograph of collapsed discharge in Neon. 
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Probe characteristics curve, Run 1, near cathode 


Fig. 12. 
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Lower plot: Intensity of detector current Ic, and 
electron density, n,vs distance, Run 2 
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Fig. 37. lon Oscillations 5 mm from 
cathode. Scale 005 Volt/cm 
vertical; 5 micro-sec/cm hori- 
zontal. Tube pressure 10 microns, 
Discharge current 65 milliamperes. 





Fig. 39. Ion Cscillations 21mm 
from cathode. Scale .02 Volt/cm 
vertical; 5 micro-sec/cm hori- 
zontal. Tube pressure 10 microns, 
Discharge current 65 milliamperes. 
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Fig. 38. Ion Oscillations 18 mm 
from cathode. Scale .02 Volt/cm 
vertical; 5 micro-sec/em hori- 
zontal. Tube pressure 10 microns, 
Discharge current 65 milliamperes. 
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Fig. 40. Ion Oscillations 33 mm 
from cathode. Scale .02 Volt/cm 
vertical; 5 micro-sec/cm hori- 
zontal. Tube pressure 10 microns, 
Discharge current 65 milliamperes. 





Fig. 41. fon Oscili@tions 28 mm 
from cathode. Scale .02 Voli/cm 
vertical: 5 micro-s_-c/em hori- 
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ion oscillation amplitude, Vi, vs distance; Run 2 
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Fig. 46 a. 





Fig. 46b,. 





Fig. 46 c. 


Fig. 46 a, b, c. Contaminated sections of the tube and internal structures 
a. Deposition of copper on glass envelope 
b. andc. Deposition of copper and copper 
oxides on tantalum anode 
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Photograph of normal glow discharge in 10 


ae | 8 microns of Neon. Illustrating the glowing - 
im Te wf region of maximum interaction. Anode to 

0 t a ee Cathode spacing 40 millimetors, Discharge 
| ils. ~~ 7% |) current 65 milliamperes probe withdrawn. 


- ~ i> Photograph taken through dildymium filter. 
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ow 





(mM) 


ANODE TO CATHODE 





1 0)4 


Dl. 


Zi 


eo 


14, 


PD « 


Lhe 


BIBLIOGRAPHY 
I Langmuir and H. Mott-Smith, Gen. Elec Rev. 449 (July-Dec 
1924) 
I. Langmuir, Phys.Rev 26, p. 585 (1925) 
Arthur F. Dittmer, Phs. Rev. 28, p. 507 (1926) 
Lewi Tonks and I. Langmuir, Phys. Rev 33, p. 195 (1929) 
H J. Merrill and Harold W. Webb, Phys. Rev. 55, p. 1191 (1939) 


E. B. Armstrong and K. G. Emeleus. Proc Inst. E. E. (London) 
96 Pr III, No. 43, p. 390 (1949) 


T. R. Neill and kK. G. Emeleus, Proc. Royal Irish Acad., 53, Sect.A 
No. 13, p. 197-219 (1951) 


R. A. Bailey and K. G. Emeleus, Proc. Royal Irish Acad.,57, Sect.A, 
Noes , pe oe. (1955) 


T. K. Allen, R.A. Bailey, K. G. Emeleus, British Journal of Applied 
Physics, Vol. 6, p. 320-322 (Sept. 1955) 


M.A. Easly, Journal of Applied Physics 22, 590-593 (1951) 


N. L. Oleson and C. G. Found, Journal of Applied Physics 20, 
416-417 (1949) 


M.J. Druyvestyn, Zeits. Phys. 29, p. 903 (1958) 


John F. Waymouth Probe Methods for Plasma Diagnostics, peer SAY 
Microfilms, Ann Arbor, Michigan 


R. H. Sloane and E. I. R. MacGregor, Phil. Mag. 18, 193, (1934) 


R. L. F. Boyd andN. D Twiddy, Physics Dept., University College 
London (5 Sept. 1958) 


K. Takayama, H. Ikegami, andS. Miyozaki: Phys. Rev. Letters 
5, 238-240 (1960) 


a2 


-— —es \ See ee | 
ee oe ee 


—_—= = Ce Ge 
+ oe — oa &-< 


a ‘ =e os + 


7 


) ome eee = 


ee 
= = 
|, ——o — — 
=m s~s 
—_ * = @ © 





ahea<£, 
Vv —, 













——_ 


ee | no 
— > -_ «= 


—- eae . * oe 


=~ lhe 
— - is a 
— | 


a 


—— << =, 





a -¢ “wee et 


— 





a7 2 


LG. 


Po. 


Zz. 


we 


Te 


ZS). 


Za. 


ZS. 


Z0.. 


ole 


28. 


W. P. Allis, Electronic Waveguides, (Microwave Research 
Institute Symposium Series, Vol. VII, Polytechnic Press, 
Brooklyn, (1958) 149-155 


D.W Mahaffey, G. C. McCullogh, A. Garascadden, 
K. G. Emeleus, Proceedings of the Royal Irish Academy, Vol. 61, 
Sect A, No. 9 (1961) 


L. H. Putnam and H. D. Collins, Plasma Oscillations in a Low 
Pressure Argon Discharge, U. S Naval Postgraduate School 
Thesis, (1961) and Phys. Rev Letters, Vol. 7, No. 3, 77-79 
(1961) 


J. D. Cobine, Gaseous Conductors, 563-564, McGraw-Hill Book 
Co., Inc., New York and London (1941) 


R. B. Brode, Reviews of Modern Physics 5, 257 (1933) 


A. von Engel, Ionized Gases, 33, Clarendon Press, Oxford 
(1955) 


D. Bohm and E. P. Gross, Phys. Rev. 75, 1851-1865, Phys. Rev. 
79 992 (1950) 


D. Looney and S. Brown, The Excitation of Plasma Oscillations, 
Technical Report No. 273, Phys. Rev. Vol. 93, No. 5, 965-969 
(1954) 


R. A. Demirkhanov, A. K. Gevorkov and A. P. Popov, The 
Interaction of a Beam of Charged Particles with a Plasma, 

Proc. Fourth International Conference on Ionization Phenomena 

in Gases, Uppsala 1959, North Holland Pub. Co., Amsterdam (1960) 


G. D. Boyd, L. M. Field, andR. W. Gould, Phys. Rev’ 109, 1393 
(1958) 


L. Alexoff, and R. V. Neidigh, Proceedings Fifth International 
Conference on Ionization Phenomena in Gases Munich 1961, 
North Holland Pub. Co., Amsterdam (1961) 


Laurence S. Hall, Probes and Magnetic Pumping in Plasma. 
(Doctoral Thesis) University of California Lawrence Radiation 
Laboratory; Published by Dept. of Commerce, Washington, D.C. 
(1961) 





29% 


30. 


Lewi Tonks and I. Langmuir, Phys. Rev. 33, p. 202, (1929) 


Robert W. Young, Frequencies of Simple Vibrators and Musical 
Scales, p. 3-102, American Institute of Physics Handbook, 


Maple Press, (1957) 





APPENDIX 

Late in this work an effort was made to establish whether electron 
oscillations could be detected near the cathode in an Argon discharge. A 
discharge was established for this purpose at a pressure of 8.7 + 0.2 
microns with a current of 90 milliamperes. 

No oscillations were detected within 5mm of the cathode at these 
values of current and pressure. Oscillations were detected near the 
center of the discharge as found by Putnam and Collins in last year's work 
These oscillations reached a maximum detector current intensity of 2 x ie 
amperes at a position 15mm from the cathode. 

We have also been able to calibrate the pressures we have reported, 
which were based on our modification of the Veeco vacuum gauge. This 
modified vacuum gauge circuit described in this paper was calibrated with 
a McLeod gauge. It was determined that the response of the Veeco gauge 
was not linear as assumed. All pressures reported in this paper as 10 mi- 
crons of Hg. must therefore be corrected to read 58 + 1 micron It was 
also determined that this non-linear response occurred for the gases Argon 
and Helium, although the factors for correction were not the same for the 
various gases. The pressure for Argon, reported above, we believe is 
accurate within the limits specified. 

This increased pressure of Neon in our tube would result in a mean 
free path for electrons in Neon on the order of 1.1 to 1.5 cm, however, the 
phenomena described in section three of the work we believe is still valid 


Further experiment in anode cathode spacing down to slightly less than 7mm 
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indicated little if any effect on the discharge characteristics. The only 
noticeable effect was a visible shortening of the plasma column which 
normally fills the greater part of the length of our tube. For these very 
small spacings the plasma surrounded the anode and cathode in the rough 


form of a sphere with a diameter of approximately 5cm. 
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